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Abstract 



We present our results of a numerical investigation of the behaviour of a system of two solitons 
in the (2+1) dimensional CP 1 model on a torus. Defined by the elliptic function of Weierstrass, 
and working in the Skyrme version of the model, the soliton lumps exhibit splitting, scattering at 
right angles and motion reversal in the various configurations considered. The work is restricted 
to systems with no initial velocity. 

1 Introduction 

Physics in (2+1) dimensions is an area of much active research, covering topics that include 
Heisenberg ferromagnets, the quantum Hall effect, superconductivity, nematic crystals, topological 
fluids, vortices and solitary waves Most of these systems are non-linear. In their mathematical 
description the well-known family of sigma models plays a starring role. One of the simplest models 
in (2+1) dimensions which is both Lorentz covariant and which possesses soliton solutions is the 
CP 1 or 0(3) sigma model. Such solutions are realisations of harmonic maps, by itself a long- 
established area of research in pure mathematics. However, analytical CP solutions have only 
been found in the static (2+0) case; their dynamics is studied using numerical methods and/or 
other approximation techniques. 

Sigma models are also useful as low dimensional analogues of important field theories in higher 
dimensions. In fact, the sigma CP model in two dimensional space exhibits conformal invariance, 
spontaneous symmetry breaking, asymptotic freedom and topological solitons, properties that re- 
semble some of those present in a number of forefront field theories in (3+1) dimensions. Amongst 
the latter we have the Skyrme model of nuclear physics Q| . Initially proposed as a theory of strong 
interactions between hadrons, it can now be regarded as a low energy limit of quantum chromo- 
dynamics The Skyrme scheme assumes that its topological solutions (skyrmions) correspond, 
at a classical level, to ground states of light nuclei with the topological charge (Brouwer degree) 
representing the baryon number. Of course, to compare with the properties of real (physical) 
nuclei one has to add to the classical results various quantum corrections. 

Planar analogues of the Skyrme model involve the addition of some extra terms to the original 
CP 1 lagrangian in order to stabilise the field solutions -the 'baby skyrmions'. Without them, 
the invariance of the pure, planar CP theory under dilation transformations would lead to the 
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instability of its soliton-lumps. In the traditional approach, where the solitons are harmonic maps 
K2 U {00} w S2 i— > S2, one adds two terms: A Skyrme-like term which controls the shrinking of 
the lumps and a potential-like term which controls their expansion. Properly implemented, this 
procedure yields stable solitons as confirmed by numerical simulations 0]. 

Lately, attention has also being paid to the CP 1 model on a torus T2 where the solitons are 
maps T 2 >— ► £2; this approach amounts to imposing periodic boundary conditions on the system. 
A characteristics of this model is that there are no solitons of topological charge one, a feature 
arising because genus (Ta)=l. 

Recent investigations have unveiled a rich diversity of phenomena in the toroidal model that 
goes beyond the two-lump and annular structures one might expect by analogy with the model 
on [|, |(| This article continues our earlier studies of CP 1 skyrmions that were found to 
undergo splitting when defined through the elliptic function of Weiestrass . 

In the following section we define the toroidal CP 1 model. The numerical procedure is ex- 
plained in section |3. Section ^ reviews previous findings and presents our new results. The paper 
closes with section % with some concluding remarks and suggestions for future research. 



2 The CP 1 model on a torus 

Our CP 1 skyrmion model is defined by the lagrangian density 
C 



\W t \ 2 - \W X \ 2 - 2\W y \ 2 
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(1 + \W\ 2 ) 2 
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9fl (W X Wy-W X Wy) 2 

+ 261 (1 + |W|2)4 ' ^ £T2 > C 1 ) 

which is the original CP 1 lagrangian modified by the addition of a Skyrme, #i-term (#1 £ K). 
We adopt the notation Wt = dtW,W x = d x W,W y = d y W for the derivatives; the bar denotes 
complex conjugation. 

In order to obtain stable lumps on the torus it is sufficient to supplement the pure CP 1 
lagrangian with a 81 term as shown in equation (jf|). This is also the case in the (3+1) dimensional 
Skyrme theory. As remarked in the introduction, in order to stabilise the lumps in the traditional 
CP 1 model on 3? 2 we have to add, other than the 81 term, also a potential-like term. 

Also it is worth mentioning that on T 2 one no longer has the problem confroted in the extended 
plane, whose non-compactness brings about formal difficulties in defining the metric on the moduli 
space of static soliton solutions Q. 

For the CP 1 model to be defined on a torus we require the complex field W to obey the 
periodic boundary conditions 

W[z + (m + in)L]=W(z), Vt, (2) 

where m,n = 0, 1, 2, ... and L is the size of the torus. The static solitons (instantons) are elliptic 
functions which may be written as 

W = Xp(z-a) + b, X,a,beZ, (3) 

where p denotes the elliptic function of Weierstrass. The partial fraction representation of p reads 

00 

p(u) = u~ 2 - ^{[w - (m + in)L]~ 2 - [(m + in)L]~ 2 }, (4) 
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the summation being over the integers m, n excluding the combination m = n = 0. A compren- 
hensive treatment of elliptic functions can be found in || B|. 

The function (Q) is of the second order, hence (^[) represents solitons of topological charge 2. 
A particularity of our instantons is that they have no analytic representative of charge one (the 
model on 3?2 has representatives in all topological classes). It is important to note that (H) is 
an approximate solution of the model ([j]), and it becomes an exact static solution in the 0(3) 
limit (#i=0) where it satisfies the ensuing field equation. This means that the CP 1 lumps should 
evolve only for 6\ ^ 0. 

Our first investigations || of periodic solitons involved the use of the Weierstrass pseudo- 
elliptic function a. Then W was taken to be given by: 

"-»n3H& x>=x>> ») 

3=1 3=1 3=1 

where the accompanying selection rule between the zeroes (c,) and poles (dj) guarantees that W 
is elliptic. 

Observe that when « = 1 in (^) the constraint between the zeroes and the poles must be 
relaxed lest a trivial configuration W is desired. Although this procedure renders W pseudo- 
periodic, starting with such a field a periodic ansatz of topological carge one was constructed in 
||. This gave a field which was an approximate solution of the equations of motion for 9\ 7^ 0, 
but became singular as 9\ — > 0. 

The power series for a on a square torus may be cast into the form 

oo 
3=0 

With the help of equation (^) the coefficients Gj can be calculated by expanding 

cr ( M ) = / [p(«) - l/« 2 ]- 



In general, the coefficients are written in terms of the so-called invariants gi (L) and gs(L). 
However, square boundary conditions effectively set gg = 0, the lemniscate case, leading to a 
relatively simple expression of the form (^). 

Let us also add that we can re-express the p function through a via the formula 

p(u) = -^]n[<7(«)], (7) 

and then perform the computation of a with the same numerical subroutine which was used in 
reference M. 

Note that each factor a(u)/a(v) in the field (^) can be used to represent a single soliton, 
providing a setting to studying more or less independent lumps in all topological classes. Solitons 
in the <r-picture may, for example, be boosted independently. On the other hand, through equation 
([]) we can only have solitons of even topological index because p is itself an elliptic function of 
order two. Lumps in the p-picture are less independent than their siblings in the cr-formulation, 
for changing the values of the parameters in always affects both lumps. Strictly speaking, 
however, truly independent solitons can only be obtained in the asymptotic regime of large lump 
separation, which really never happens on a compact manifold Jlo| . 

Now let us display some useful relations satisfied by the p function, relations which follow 
from (0). They are: 

p(-u) = p{u), p{iu) = -p(u), p(u) = p{u). (8) 
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Moreover, there exists a useful algebraic relation between p and dp/du on a square torus which 
reads 

[^] 2 = Ap{u)[p{uf - p(L/2)]. (9) 

It can be deduced from the equations ([]) and ([]) that the function p is purely imaginary on 
the diagonals bisecting the fundamental cell and real on the central cross and on the boundary of 
the cell: 

{imaginary : on central diagonals; 
(10) 
real : on central cross and boundary. 

The static energy density associated with the field (^) can be read-off from the lagrangian 
density (Q). Using formula (g) we may write 

Pictures of E reveal a distribution of lumps localised in space. The parameter A is related to the 
size of the lumps, b determines their mutual separation and the parameter a merely shifts the 
system as a whole on the torus. Throughout our simulations we have chosen the values 

A =(1,0), a = (2.025,2.05), 6i = 0.001 (or zero), 

and studied different configurations by varying b. 



3 Basic numerical set up 

We have taken fields of the form (^) as the initial conditions for our time evolution, studied numer- 
ically. Since the field W may become arbitrarily large, we have preferred to run our simulations 
in the ^-formulation of the model. Its field equation follows from the lagrangian density (^) with 
the help of the stereographic projection 



W + W . -W + W \W\ 2 - 1 

^]wY+i ,l WW+i \w\ 2 + r' 



(12) 



with the real scalar field <j> satisfying <j}.<j> = 1. One has 
= (a M 9 M - $.d"dj)$ 

- (dj.d„$)(d^.d"$)$], 



(13) 



where /i, v = 0, 1, 2 are the Lorentz indices; as usual we have t,x,y — x°, x 1 , x 2 . 

We compute the series (^) up to G5, the coefficients Gj being in our case negligibly small for 
3 > 5: 



Go 
Gx 

G-2 

G 3 
Gi 
G, 



= 1 

= -0.7878030 

= -0.221654845 

= 9.36193 x 10" 3 

= 7.20830 x 10~ 5 

= 2.37710 x 10 -5 



We have employed the fourth-order Runge-Kutta method and approximated the spatial deriva- 
tives by finite differences. The laplacian has been evaluated using the standard nine-point formula 



4 



and, to further check our results, a 13-point recipe has also been used. Respectively, the laplacians 
are: 



V 2 = 



14 1 

4 -20 4 
14 1 



6 x a 2 





-1 
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12 


1 


12 


-48 


12 


1 


12 


1 




-1 




10 x a 2 



The discrete model has been evolved on a n x x n y — 200 x 200 periodic lattice with spatial 
and time steps 5x—Sy— 0.02 and <5t=0.005, respectively. The vertices of the fundamental lattice 
we have used for our simulations were at 

(0,0), (0,L), (L,L), (L,0), L = n x xSx = 4:. (14) 

Unavoidable round-off errors have gradually shifted the fields away from the constraint <f>.(j> — 1. 
So we rescale 




every few iterations. Each time, just before the rescaling operation, we evaluate the quantity 
/i = <f).(j) — 1 at each lattice point. Treating the maximum of the absolute value of /x as a measure 
of the numerical errors, we find that max|/x| » 10 -8 . This magnitude is useful as a guide to 
determine how reliable a given numerical result is. Usage of an unsound numerical procedure in 
the Runge-Kutta evolution shows itself as a rapid growth of max|/it|; this also occurs, for instance, 
in the limit 6\ — > when the unstable energy lumps become infinitely spiky. 



4 Results 



In reference |J we considered two cases: b=(0,0) and fe=(l,0). For 6=(0,0) one has a most evenly 
spread energy distribution. Notably, this homotopy-two class confi gur ation exhibits four lumps 
(rather than two) sitting on the central diagonals of the basic grid (|14[). (See the top-left graph 
of figure [j].) When evolved from rest in the Skyrme scheme, this quartet has moved along the 
diagonals under the action of a net repulsive force. They have evolved in a manner that resembled 
a scattering at right angles as illustrated in upper half of figure ^. No splitting has been observed. 



The state 6=(1,0) corresponds to a couple of lumps placed along the cordinate axes as seen in 
figure (top-right). This is certainly a more familiar picture for solitons belonging to a charge-two 
topological sector. Unfamiliar, though, is their novel dynamics: starting from rest after a while 
each soliton splits into two lumps. The split is in the direction perpendicular to the line joining 
the solitons. As time goes by, the offspring skyrmions glue back together, split again and so forth. 
This is depicted in the lower part of figure ti and is the splitting phenomenon described in the 
introduction. The time at which the lumps begin to split is t ~ 7. 

The lower half of figure |l| shows plots of E max vs. t for the above cases, including the unstable 
0(3) case. In the latter, as expected, the lumps do not move at all (let alone split) with the 
passing of time -as long as the initial speed is zero. This is in accordance with our expectations, 
for, as we recall, the field (Q) is a static solution of the pure 0(3) model. 

Note that in general b — (a, 0) corresponds to solitons initially located on the central cross of 
the grid. If a > (a < 0) the lumps lie on the vertical (horizontal) axis. Our qualitative results 
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are unaffected by reasonable values of a. In connection with these configurations it is interesting 
to consider the case where the simulations run in the Skyrme format up to t = to, after which 
they evolve with d\=0. In other words, we perform the simulations with a Skyrme field as the 
initial condition for the 0(3) evolution. Since the splitting forces act only when 9i ^ 0, we expect 
that the larger to is the sooner the solitons will begin to divide up. 



Let us analise the situation arising from b = (— 1, 0) , which positions the lumps along the 
horizontal axis. Figure |^ illustrates the evolution of the system for the sample cases to= 1-3, 1.5. 
We find that after splitting the extended structures reunite but eventually break up for good. In 
agreement with our estimates, for to=1.3 the system splits later than it does for to=1.5. Also 
apparent is that the lumps show a stronger tendency to glue back in for to=1.3. We may now 
compare these plots with figure ^, which exhibits the splitting situation for b — ( — 1,0) in the 
limit to = oo, i.e., when the routine evolves with the Skyrme term on at all times. 

The previous diagrams suggest that a to might exist for which the lumps would come back 
together without further division. After some trial and error we have found that such critical 
time is approximately to w 1.25. The plots presented in figure ^ show that the forces brought 
about by the skyrmionic initial conditions are just enough to set the solitons oscillating in a 
break-up-join-up fashion. The skyrmions eventually settle together and stay that way for as long 
as the numerical procedure can be trusted: In the pure format E max shrinks indefinitely as usual 
although it takes quite a while for the system to blow up [t ~ 100). This soliton is almost a static 
solution of the model. 

Let us point out that when defined via equation (^) the periodic chunks of energy were found 
not move at all when vo = 0, neither in the pure case nor in the Skyrme case ||. Such result is 
rather unexpected for skyrmions, for they are only approximate solutions of the field equation. 
Therefore, solitons obeying (^) look sensitive to the choice of function describing them. Not so on 
K2, where the qualitative behaviour of the solitons is qualitatively the same for all the choices. 



Let us now study 



fe=(0,l) . For this value the energy distribution (^lj) associated with (g) 
has the form of two lumps placed on a bisecting diagonal of the elementary grid (figure ^). Our 
numerical simulations show that the skyrmions attract each other and collide at the centre of the 
grid. They coalesce indistinguishably for a moment and then re-emerge perpendicularly to the 
initial line of approach, scattering off at 90° . Afterwards they continue towards the corners (which 
are the same point in a periodic set-up like ours) where they again scatter off at 90° . These events 
are depicted in figures [7 and [| the corresponding E max (t) diagram is the Skyrme curve shown in 
the lower half of figure 9j. Such multi-scattering process goes on indefinitely, and no splitting is 
observed. 

We have also evolved the b = (0, 1) case with 61 — 0. This time the energy chunks have 
not moved at all with the passing of time until they have become too narrow and the numerical 
procedure has broken down (see figure |^) . The shrinking of solitons and the problem of singularity 
formation in this model was predicted in using the geodesic approximation. 



Next we consider 



6= (1,1) . Here we have a pair of lumps initially situated at some almost 
random points, neither on the central cross nor on the central diagonals of the lattice. Figures [Hi] 
and [ll] reveal an interesting evolution for these skyrmions. As time elapses they move along the 
flat torus along the x-axis, disappearing/re-appearing through the lateral edges of the basic cell 
( |l4| ) in a continuos, periodic motion. For instance, the skyrmion at the upper (lower) half of the 
grid journeys towards decreasing (increasing) x whilst keeping its y coordinate roughly constant. 
The lump disappears into the line x — (ie=4) re-emerging from the opposite side x = 4 (x = 0). 
By imagining the flat torus as the product manifold of two circles, T2 = Si X Si, we can visualise 
the motion just described as the trajectory on a circle with our x coodinate corresponding to an 
angle. 

The two lumps continue moving along and, around t=69.5, they reverse their motion: The 
skyrmion at the upper (lower) half of the grid starts towards increasing (decreasing) x whilst 
keeping its y coordinate roughly constant. Now the upper (lower) lump will disappear through 
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the line x — 4 (x=0) and re-emerge on the opposite side x = (x = 4). Figure |l2| presents the 
variation of the coordinates in time, as well as the evolution of the peak of the energy density. No 
splitting was observed for fields of the type 6=(1,1). 

Thus, for skyrmions whose initial positions are determined by b = (a, 0) our results indicate 
that 

• splitting only occurs for b = (a, 0), that is, when the solitons originally lie on the 
central cross of the network (|14|). A repulsive, splitting force acts within each lump. By 
using a Skyrme field as an initial condition for the 0(3) evolution, it is possible to introduce 
some attractive forces so that the lumps stick back together. 

• Head-on right angle scattering only takes place for b — (0,/3), that is, when the 
skyrmions initially sit on the central diagonals of the grid. There is a net attractive force 
between the lumps. Now, 

— if j3 7^ then we have two lumps located on a diagonal; 

— if P = then we have four lumps symmetrically situated on both main diagonals. In 
this case, the scattering becomes a breather-like vibration with the four lumps moving 
in towards the centre of the grid and moving out again. Interestingly, a look at the 
property ( p^| ) indicates that for fe=(0,0) the soliton W = p{z ~ a) is purely imaginary. 

• Neither splitting nor scattering is observed when b = (a,0), a, (3 ^ 0. Here we have a 
configuration of two energy humps propelled along the lines parallel to the edges of the cell 
and, later on, experiencing a 180° change in their sense of motion. This looks like scattering 
at a large impact parameter. We remark that all systems considered in the present paper 
start off from rest. 



The right angle scattering is a generic feature of sigma-type models [n2|. In the case of the 
CP 1 dynamics on T2 the scattering at 90° was theoretically hinted at in |7[ (pure case), where 
the initial value problem was defined in terms of the p function (Q); the solitons were evolved 
using the geodesic approximation. In reference jO] we evolved the p solitons using a numerical 
simulation of the full CP 1 model (pure case) and confirmed the 90° scattering predictions of 
Earlier on ||, we had observed the said scattering on T2 (for both the pure and Skyrme cases) 
with solitons expressed through the a function <Mj. The present work further establishes these 
results, with the added interesting feature that the lumps, in the Skyrme case, collide and scatter 
off even though they have no initial speed. 

Among other sigma models that exhibit 90° scattering we can mention monopoles and 
vortices JL4| . 

Finally, it is worth to note from the bottom-left diagram of figure |l| that the lumps for the 
b = (0, 0) case shrink very little, which suggests quite a stable configuration. So we have performed 
a number of simulations with b—(0,0) in the pure 0(3) format. Indeed, as figure jl3| shows, the 
solitons shrink at a very slow rate and can be regarded as stable for practical purposes. Contrary 
to the case of Si =0 processes with no initial speed, a displacement in this four- lump configuration 
is observed, one of breather-like characteristics. 
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Figure 1: Energy density configurations at t — and the evolution of the maximum E max . Above 
we see the special case featuring four lumps [6=(0,0)] and the more common situation with a pair of 
lumps [b=(l,0)]. The evolution of E max is illustrated bellow, both for the Skyrme and pure 0(3) case. 
The instability of the latter is manifest by the lumps shrinking non-stoppingly. 
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Case b=(U,U) b density tor t=3U 




12 3 4 

Figure 2: Above: Skyrmion for the exceptionally symmetric case 6=(0,0). The four lumps stay on 
the diagonals and move towards the corners and coalesce. They break-off again and proceed back to 
the centre of the lattice at £=51. Below: The b=(l,0) skyrmions split each in two lumps that transit 
complicated paths; the labels a — h refer to one of the 'half-lumps'. The t—30 picture features the 
situation shortly after the 'fractional' lumps reunite at d (and at its symmetrical point) and begin to 
travel centrewards. 
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Figure 3: The CP 1 solitons b = (— 1, 0) are run with B\ ^ (Skyrme format) up to t = to, at which 
time 9\ is set to zero (pure format). As expected, the lumps divide up more readily for larger t$. The 
change in the vertical coordinate of E m ax as time goes by, and the trajectory plots for the two cases 
to = 1.3, 1.5 are shown. 
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11 




12 




Figure 6: Energy distribution, equation ([il|), for 6=(0,1). The extended structures sit along the main 
diagonals of the lattice. Starting off from rest, these entities attract each other, collide at the centre 
of the grid and scatter at right angles with respect to the initial direction of motion. (See figured.) 
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t=6,7 p 19.5, 22 




Figure 7: The solitons of picture || collide at the centre (i=6) and scatter off at right angles (t =7) 
A similar event occurs at the corners (t = 18.5), after which the solitons proceed back to the centre 
Such scatterings take place time and again, in periodic cycles. No splitting is detected. 
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12 3 4 



Figure 8: Trajectory plot further illustrating the 6=(0,1) lumps described in the previous picture, 
figure [t]. The Skyrme term introduces forces that make the lumps cruise along the diagonal and collide 
head-on (t=6). The scattering at ninety degrees is apparent: Around <=18.5 the skyrmions bump into 
each other at the corners (x, y)=(0,4),(4,0), re-emerging perpendicularly through (x, y)=(0,Q),(4,4). 
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Emax=448.7 at t=1 0.5 



The lumps breakdown (t > 10.5) 




Figure 9: Pure CP 1 case corresponding to the value 6=(0,1). The lumps stay still at their initial 
positions (top-left) as time elapses. They break down as illustrated at the top-right diagram. Below 
we have the evolution of E max , including the Skyrme case. 
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Emax=1 1.2127; t=G 



Emax=15J0951 >23 



Emax=1 23873; t= 41 

E7 



Figure 10: Showing the dynamics arising from b = (1,1). Although no initial speed is impinged on 
this system, the skyrmions move resembling a collision at a certain impact parameter. One skyrmion 
moves towards decreasing x whereas the other skyrmion proceeds in the sense of increasing x. The 
disappearance/appearance of the lumps through the edges of the lattice reflects the periodicity of the 
network (111). 
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Figure 11: Motion reversal: The skyrmionic structures of figure |lOj continue their itinerary but, at 
t = 69.5, their motion begins to reverse. 
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10 20 30 40 

t 

Figure 12: Above: The evolution of the coordinates of E max corresponding to the event depicted 
figures |l^ and [ll} Below: The evolution of E max itself. 
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5 Concluding remarks 



The CP 1 model in (2+1) dimensions has a very rich structure. On the torus we have uncovered 
several characteristics qualitatively different from those seen in the usual model on the compact- 
ified plane. For example in the two skyrmion sector there exits a configuration featuring four 
lumps rather than two. (Clearly, care should now be exercised when referring to the topological 
charge as the 'lump number'.) We also have the interesting fact that no one soliton solutions exist 
on Ti. The lumps are not stable in the pure 0(3) set-up, but become stable in the Skyrme model, 
constructed by the addition of just one extra term to the 0(3) lagrangian. 

Skyrmions on T2 have many other important properties. For instance, their behaviour depends 
on the elliptic function used to define them. Thus, skyrmion fields expressed in terms of the a 
fields, e.g., equation (jE|), evolve differently from those expressed through p, e.g., equation (^). In 
the former case, energy chunks started off from rest stay still in their initial positions as time goes 
by. In the latter case, we have encountered novel phenomena like lump-splitting, scattering at 
ninety degrees and motion reversal, despite no initial speed being given to the systems. One of 
the aims of this research has been to investigate various properties of solitons on T2 by employing 
alternative elliptic functions. Comparison of the various solitons generated by alternative elliptic 
functions (and with the lumps on K2) may help us to gain a deeper insight into the CP 1 dynamics. 

One of our main results is the existence of a mode of splitting a skyrmion into two lumps 
(though we have to work with two skyrmions which then split into four lumps). There clearly is 
a potential barrier which, when overcome, allows the skyrmions to divide up. We have managed 
to overcome this barrier by starting our simulations with a Skyrme initial condition which runs 
up to a t = to, when the Skyrme term is switched off: The evolution continues with 6\ — 0. Such 
an approach adds a little energy to the system and then transforms it to the mode responsible for 
overcoming the barrier. This barrier is a numerical artifact, brought about by the discretisation 
procedure. 

The longer the system runs with 6\ 7^ the more energy is transferred; hence there seems to 
be a minimal time of such a simulation t — to below which the skyrmions do not break up. This 
we have seen in our simulations - our estimates of to gave us to ~ 1.25. 

In the present paper we have also learned that splitting takes place only when the parameter 
b in W — p + b satisfies b £ K — {0}. For some values of b we have found that the forces operating 
in p-lump systems with vo — either lead to head-on collisions and a subsequent 90° scattering 
or lead to a motion reversal. This property is qualitatively different to any other CP 1 structures 
that we know of. 

The table below summarises some of our results: 



Im(b)/Re(b) 


= 


7^0 


=0 


4 lumps on diagonals; 
breather 


2 lumps on central cross; 
splitting 


^0 


2 lumps on diagonals; 
90° scattering 


2 lumps elsewhere 
motion reversal 



In future work we hope to report on: 

• configurations with vo 7^ on the central cross to further analise the splitting phenomenon; 

• lumps initially situated slightly off the central cross [b = (a, /3), < f3 <C 1]; 

• lumps slightly off the main diagonals [b — (a, P), < a < 1]; 

• obtain deeper insight into motion reversal. 

A look a bit farther ahead will necessarily include configurations with higher Brouwer degree, 
as well as states defined through other functions, e.g., the elliptic ones of Jacobi and the pseudo- 
elliptic theta functions. 
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